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Abstract 

Neutrino oscillations induced by a flavor-dependent violation of the Einstein Equiva- 
lence Principle (VEP) have been recently considered as a suitable explanation of the solar 
i^e deficiency. Unlike the MSW oscillation mechanism, the VEP mechanism is dependent 
on a coupling to the local background gravitational potential $. We investigate the dif- 
ferences which arise by considering three-flavor VEP neutrinos oscillating against fixed 
background potentials, and against the radially-dependent solar potential. This can help 
determine the sensitivity of the gravitationally-induced oscillations to both constancy and 
size (order of magnitude) of "I". In particular, we consider the potential of the local super- 
cluster, \4>sc\ = 3x 10~^, in light of recent work suggesting that the varying solar potential 
has no effect on the oscillations. The possibility for arbitrarily large background potentials 
in different cosmologies is discussed, and the effects of one such potential ($ — 10~^) are 
considered. 



1 Introduction 

A resurgence of flavor-oscillation solutions to the Solar Neutrino Problem 
has recently inundated the world of particle astrophysics. Of these, the 
two main competitors are the tried and true MSW mechanism ^, and 
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the more radical VEP (for "Violation of the Equivalence Principle") oscil- 
lation model Ip. The MSW mechanism is often the more seriously consid- 
ered one, due to the basic underlying hypothesis of VEP, which requires a 
generation-dependent violation of the Einstein Equivalence Principle (i.e. 
flavor-dependent coupling to the background gravitational potential in 
the neutrino sector, and tends to give most physicists a headache. Con- 
versely, there is no experimental justification for any flavor of neutrino to 
have mass, so MSW can also be considered a large assumption. 

Another consideration is the number of flavors with which to work. Many 
studies of neutrino oscillations ^ take the simpler two-neutrino approach 
by only considering the effects of a flavor eigenstate \i^w) = {^e, I'/i)'^ ■ There 
are at least three flavors known in nature, though, and various studies g, 0] 
have suggested that the two-flavor model short-changes the viable solution 
set by considering only a minute portion of the overall parameter space. The 
authors of |^ suggest that three-flavors are sufficient to solve the SNP. 

Proper experimental results from solar neutrino detectors could help 
isolate the individual values of the mass eigenvalues and mixing angles which 
contribute to oscillations. Experimental detection of oscillations induced 
by the MSW mechanism can yield determination of the neutrino energy 
eigenvalue differences Amfj/{2E), and hence the neutrino masses rrij. If 
VEP is the true underlying mechanism, the analogue of the energy eigenvalue 
difference is 2E^Afij. As long as we are aware of the background potential 
$ that causes the oscillations, we can isolate the values Afij, and thus the 
individual fiS. While experimental verification of MSW can teach us new 
neutrino physics, a detection of VEP oscillations can additionally yield new 
information on solar astrophysics and General Relativity (the universality 
of the equivalence principle) . 

For solar neutrinos, one would generally expect $ = ^©(r), i.e. that of 
the solar interior. However, there is some debate as to whether or not the 
propagating neutrinos feel the effects of a potential difference, or whether 
they feel the overall effects of the strongest potential [|9|. If the latter is 
the case, then the gravitational potential of the local supercluster {(pscl = 
3 X 10~^ overwhelms the former by almost a factor of 10. In effect, the solar 
neutrino problem has actually become the solar neutrino problem: without 
proper knowledge of the background potential which contributes to the flavor 
oscillations, we can only learn as much new physics as we could in the MSW 
case. That is to say, we can only determine that neutrinos do oscillate, but 
not their gravitational eigenvalues. 

This then raises the question of how much of a difference the choice of 
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potential makes in the overall process. The following paper will examine the 
dependence on potential energy of flavor oscillations in a three-generation 
framework, from two specific standpoints. To begin with, we will compare 
the effects of the radially-dependent solar potential (pQ{r) with a constant 
one of equal magnitude (the solar surface potential, \(I)q{Rq)\ = 2.12 x 10~^ 
[jiol ) . This will give an indication of the sensitivity of the data to variations 
over a fixed magnitude. Secondly, a comparison will be made of several 
constant potentials, to examine dependence on varying magnitudes. These 
will primarily include the solar surface potential, and the local supercluster 
{(/^scl = 3 X 10^^. An arbitrarily large potential of the order $ = 10~^ will 
be considered, as well, for reasons discussed in section H 



2 MSW and VEP Oscillation Mechanisms 

At heart, the MSW and VEP mechanisms are the same. Both rely on the 
existence of two distinct eigenstates of the neutrinos. In the former case, 
there are mass and electroweak eigenstates, \i')m and \i^)w- For the latter, 
the mass eigenstate is replaced with a gravitational one, \i')gi as the VEP 
neutrinos are considered masslessQ. The equations of motion are diagonal 
in the mass/gravitational (M,G) basis, 

i-^HM,G = H\v)m,g , (1) 
ar 

where @ is derivable from the Pauli equation [Q]. To investigate neu- 
trino flavor dynamics, the (M/G) states \i')m,g are rotated to the elec- 
troweak basis via |z^)vi/ = ^|i^)m,g j ^ ^ SU{N) for flavor^. Due 
to charged current interactions in the sector, an additional interaction 
A = diag{\/2GFNe,0,0, ■■■) must be added, giving the final equation of 
motion as 

^-^W)w = {VHV-^ + a] \u)w = H'\u)w ■ (2) 

for neutrinos of energy E. Here, A{r) = diac/(2\/2GFA'e(r), 0, 0, ...), with 
Fermi's constant Gp, and the radially-dependent solar electron density 
Ne{r). Subtracting an overall factor of (and hence an unobservable 

^ There has been some investigation of the possibihty of massive VEP neutrinos, but 
those are not considered here. See for further information. 

^Although V is a four parameter matrix (3 real rotations and a complex phase), we do 
not consider CP violations here, and thus eliminate the complex phase. 
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phase in the wavefunction) from (|2|), and re-diagonahzing, one obtains eigen- 
states dependent of either values of Am?j/2£' or 2E\^\Afij, i,j / l,i 7^ j. 
The main signature of VEP is its inverse dependence on the neutrino energy 
E, as compared to MSW. 



3 Three— Flavor VEP Formalism 

A brief review of the underlying formulation of the VEP mechanism is re- 
quired before proceeding. Three-flavor oscillations take place due to the 
existence of two non-trivial spinor eigenbases, = (^^i, 1^2, 1^3)"^ > W)w = 
(fe, Vfj,, i^t)'^, related by a rotation matrix V3 G S0{3) {i.e. no CP violations), 

W)w = V3W)g ■ (3) 

The gravitational eigenbasis evolves according to the massless Dirac 
equation 

ipi^G = , tl;G = e^p{-iEkt} , (4) 

k 

The Hamiltonian derived from Eq.(Q) becomes off-diagonal under the 
change of basis (|3|) , 

H' = V3H{<!>)Vi^ +A{r) , (5) 

with 



H{<^) = 2E\<P\ 



/ /i 
/2 



V /3 

A{r) = 2V2GFNeir)diag{l,0,0} . (6) 

We explicitly express the dependence on the potential of H = H{^) 
for our later consideration of the various forms of <I> = {(j)Q{r),(f)sc,4'SD}, 
but omit it for the other Hamiltonians (and assume its dependence to be 
understood in these cases). The matrix A{r) is introduced by the charged- 
current interactions of the i^gS (hence the triviality Aij = V {i,j} 7^ {1, 1}), 
and results in a shift in the energy eigenvalues E^ from their values in ^. 
A new basis \i')mg and set of corrected eigenvalues {El,E2,E'^} may be 
found by rediagonalization 
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H" = V^H'{V^)-^ , (7) 

where V3" = {612, (^1^,^23) corrected rotation matrix, with 

matter-enhanced angles (see [R] for the definition of these, or |12, 13| for 



their MSW equivalent). 

The wave equation may now be solved, subject to the modified Hamilto- 
nian H", and the subsequent survival probabilities for neutrino eigenstates 
can be found. For solar neutrinos, the value of interest is {P{i^e ^e))i 
since most detectors to date are sensitive only to first-generation neutrinos 
(this is not true of future detectors such as Superkamiokande or SNO, which 
can detect neutral current interactions from all flavors). Averaged over the 
Earth-Sun distance, the expression obtained is 



3 

2 



(P(^e-^e)) = E l(^3)l,.n(PLz).,f K^s'")!, 



Cml2Cml3 {(^ -^'1)^12^13 + PlSi2c\'i^ 

+4l2cL3 {A(1 - P2)cl2cl^ + (1 - Pl)(l - P2)4c?3 + ^2^} 
5^3 {PlP2C?2C?3 + ^^2(1 " A)s?24 + (1 " ^2)^} • 



+ 



which is a function of the matter-enhanced mixing angles , 9^ and 
eigenvalues differences Ai<2i,Ai<3i, and the regular (vacuum) equivalents 
^12,^13, and A/21, A/31. Note that {P{i'e — > ^e)) does not depend on inter- 
actions between the second and third flavors, f ^ — > z/,- , i/^ — > f ^ . 

The results of 0, |l|] offer insights into the fundamental theoretical and 



observational differences between the VEP and MSW mechanisms. The 
overall conclusion was found to be that the addition of the third neutrino 
flavor can significantly enhance the overall allowed parameter space from 
that of the two-flavor model. Unique experimental signatures in observed 
Ue fluxes can also arise between the VEP and MSW, which could be used to 
differentiate between mechanisms. 



(8) 
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4 A Potential Problem! 



While most papers on VEP tend to agree on the form of the oscillation 
mechanism^, there is some debate as to which background potential is at 
work. Most works have assumed that neutrinos feel the potential difference 
of the sun from their creation points, to their exit on the surface, which 
ranges between ~ (—7,-2) x 10~^ 0]. Alternatively, Q suggests that the 
neutrinos feel instead the constant background potential of the Local Super- 
cluster (Great Attractor), which in the neighborhood of the solar system is 
\4^sc\ = 3x 10^^. This overpowers the surface potential of the Sun by almost 
a factor of 10. Results from the COBE data suggest that the dominant local 
background potential is very close to this value p^ j. 

Furthermore, while the current dominant potential in the local neigh- 
borhood seems to be that of 05C) one should be careful to consider the 
overall potential contribution of the Universe. Following [16|, the local back- 
ground potential can be expressed as the sum of all contributing potentials, 
^locai = </>© + (t^odot + (t>MW + 4>sc + <pu + C . The contributors of the po- 
tentials are, respectively, the Earth, Sun, Milky Way, Local Supercluster, 
the Universe, and an arbitrary constant C. This is different from Equa- 
tion (3) of 1 16 1 in that the (dominant) potential (psc has been added. The 
constant C comes from an arbitrary cosmological model that might be at 
work, but is undetectable by our standards. That is, as expressed in p6| ], 
"we cannot 'step outside the universe' to [measure the actual value of C]"0 
Hence, a future cosmology might contain a C such that ^locai ~ C. If C is 
sufficiently large, it will overwhelm the other contributors. For this work, 
we will assume the value of C = 10~^, and denote this by 4>sD- 

From solar neutrinos alone, we can only hope to isolate the values of 
Afij^ (see, for example, and not the individual eigenvalues A/jj. Ac- 
cording to the aforementioned product, the values of and of Afij are 
inversely linked. So, by considering the effects of the Local Supercluster 
over the effects of the solar potential (pQ^r), one is effectively shifting 
the values of the violation parameters down by an order of magnitude. At 
first glance, this would seem to not be a big problem, but it should be noted 
that shifts of this magnitude can strongly affect the oscillation behavior for 

^ There is a brief mention in of an alternative form of equivalence principle violation, 
based on neutrinos coupling to V$, instead of "1>, in some string theories. 

^This can be thought of as a type of philosophical argument, along the lines of "how do 
we know the reality is not just an illusion?" , but with actual physical evidence to support 
it! 
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critical values of mixing angle and EEP-violations from the other flavor. 



5 {P{TJe ^ T^e)) Curves and Surfaces 

For two-flavor analyses of neutrino oscillations, it is useful to study the av- 
eraged survival probability curves {P{ve — > ^e))- Such plots are generally 
expressed as functions of the argument EAf, scaled to an appropriate value. 
These curves indicate the overall behavior of {P^i^e ~^ ^e)) with respect to 
values of violation parameters A/ [i.e. fi and /2) and mixing angle 9. In 
particular, it is of use to examine the small and large-mixing angle limits. 

For two flavors, the resonant boundary for solar neutrinos is determined 
by equating the resonance density 

_ V2\(t>eir)\EAfcos20G 

[lyejres — , 

Since the density of the sun is strictly decreasing from its center, then if the 
resonance density for a neutrino is greater than the density at which it was 
created, it will never undergo matter-enhanced oscillations and will propa- 
gate as in vacuum (note that they may still undergo vacuum oscillations). 
Using the radial solar potential (j)Q{r), one obtains a resonance boundary of 
A/ ~ 2 X 10~^^ for neutrinos of energies Ei, > 0.2 MeV 

An extension to a three flavor model requires something extra, though. 
Instead of just one specific curve for O12 (the equivalent of the two-flavor 
angle 0), there can be a whole family of curves {(P(fe — > t'e))} corresponding 
to the same 9i2- These are, of course, determined by different values of ^13, 
and imply that the large and small limits of ^13 must also be considered. 
Hence, it becomes apparent that an examination of a probability surface is 
in order. 

These surfaces are of use in three-flavor oscillation dynamics as the 
demonstrate the full range of behavior of {P{i^e — *■ z^e)) foi" all values of 
A/ji$ for a given set of ^12, ^13- This includes regions for which neu- 
trinos with resonance densities above the maximal solar electron density 
{N(,)lnax =WOgcm~^. Again, if the resonance density for a neutrino i ex- 
ceeds this value, the neutrino will not feel matter-enhancement, and will 
propagate as in vacuum. The resonance density for each flavor is deter- 
mined by the generalization of ^ to three flavors, by replacing A/ with 
A/jj for the i^i — > i^j resonances. 

Note that the region encompassed by these surfaces includes areas where 
A/21 > A/31. This corresponds to a broken hierarchy in the eigenstates 
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(i.e. /2 > /s). We would like to think that nature would follow a hierarchial 
ladder between generations, but there is nothing which physically rules out 
the possibility of breaking this. Indeed, the behavior of {P{i^e ^e)) does 
exhibit noticeable changes when the hierarchy is broken. In particular, this 
generates a change in the values of the matter-enhanced mixing angles ^21 
and 6y^. This change is less apparent for small vacuum mixing angles, while 
it can create sizable differences for large vacuum angles 

It should be noted that these plots show the survival probabilities for 
z^eS0, the easiest of the solar neutrino spectrum to detect. The majority of 
the events recorded in the Homestake i^'^Cl), as well as all events detected 
by Kamiokande, are of this type, due to their extremely high ranging energy 
spectrum of [0, 15] MeV. The figures herein have been calculated using the 
most recent solar neutrino data and Standard Solar Model data available 



[1|, g. 



The expression for {N^.y^'^ in doj) is clearly a function of the background 
potential (j)Q){r). Rearranging this, it is obvious that 

A/ = (^'^)"^"^^^ , (10) 
V2E\,j>Q{r)\ 

i.e. the resonance boundary varies inversely with the potential. So, there 
should be a decreasing limits on the values of A/ which permit resonance as 
the potential increases. For the cases we concern ourselves with here, these 
boundaries should approximately (within an order of magnitude) occur at 
A/ ~ 10"^^ for (l)sc and A/ ~ IQ-^^ for (I)sd- 

In addition to the u^s, the other main type of detectable neutrino is 
the low energy pp-neutrino, resulting from the proton-proton chain reaction. 
These are emitted over an energy spectrum ranging to 0.42MeV. Since their 
maximal energy is at least one order of magnitude different from the maximal 
^B neutrino energy, then the overall product EAfij can shift by the same 
order of magnitude. It was observed in Q that this shift can result in 
different resonant behavior between the reaction types of solar neutrinos, 
such that some can exist under the resonant barrier (2£'A/jj/10^^^MeV~ 
10^), while others can be over. The crux of this is that certain neutrino 
types will be subject to resonant behavior (for lEAfu < 10^, in units of 
lO^^MeV"^), while certain types will propagate as in vacuum {2EAfii > 
10^). While this analysis is restricted to ^B neutrinos, the interested reader 

^These originate from the thermonuclear decay reaction ^Be + Ve\ for more 

information on solar neutrino-producing reactions, see p^ ]. 
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should reference Q for an appropriate discussion. 

6 Discussion of Figures 

As previously mentioned, this work considers two distinct cases of survival 
probabilities, as a function of background potential energy. We first con- 
sider the dependence on (-P(fe Ve)) curves as a function of variability in 
potential over a specific order of magnitude, and secondly as a function of 
variability of order of magnitude. 

6.1 (j)Q{r) v.s. 0o(/?o) 

It was noted in that there is a strong influence from the ^13 contribution 
of the third flavor. In particular, for large ^13, the survival probabilities for 
Ue become effectively {P{ve ^e)) ~ s\v? 61^. This is clearly evident from 
Figures Q,|2[ which demonstrate the behavior of {P{ve ^e)) from solar 
neutrino data in the small and large O12 range, for small ^13. The end result 
of the difference in potential presents itself as a horizontal shift to the right 
for {P{i'e ^e)){4'&{P&)) when the slices are considered. The shape of the 
curve is preserved, which suggests that the change in potential has merely 
shifted the resonance density by an overall constant. Both curves assume 
the same values before and after this barrier. The portions of the curve 
for which 2EAf2i > 5 x 10^ correspond to broken hierarchy between the 
violation parameters A/ji,i = 2,3. 

Meanwhile, the story for small ^13 is set to a somewhat different tune. 
We observe vastly different behavior for both curves in each case of small 
and large 6*12 (Figures ^,^) . Both plots show that the the survival probability 
(P(fe — ^ ^e))i4'QiP-Q)) is greater than {P{i^e ^e)){<pQ{^))- The shift in 
the location of the resonance boundary 

6.2 Order of Magnitude Dependence 

In contrast to the previous section, there is a noticeable difference in (-P(z^e — *■ 
Vg)) curves and surfaces when the order of magnitude of the effective poten- 
tial Phi is considered. This shift has large implications in particular when 
it comes to deciding which potential is at work, if the VEP mechanism is 
correct. Indeed, the evidence presented here could help point to the cor- 
rect choice, subject to experimental verification {e.g. flux spectra can be 
directly computed with knowledge of Psurv). 
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For the purposes of this discussion, both the survival probabiUty curves 
and surfaces wih be examined for select choices of the potential. As in 
Section |6.1 | , the small and large 9i2 (fixed 6*13 large) curves are presented 
in Figures |5|,^. The chosen potentials here are 4>Q{r), 4>sc, and (psD- It is 
assumed that (j)Q{RQ) will vary as discussed in the previous section, and 
thus is not considered here. 

Figure ^ represents the small-angle mixings for both 12- and 13-resonances. 
These curves correspond to a choice of A/31 = 10~^^, at a fixed value of 
2£'A/3i = 2 X 10^. In this case, each potential induces a drastically different 
shape for each curve. The choice of 2EAf3i places the Solar curve in the 
double-resonance region (the can resonate to both i/^ and i^r), while the 
(j)SD curve allows only — > resonances. In fact, the shape of this curve 
approaches that for two-flavor non-adiabatic resonance [0]. By referencing 
Figure 0, it is evident that the plots for each $ represent equivalent slices 
for 00 (r) at different values of 2£^A/3i (in this case, roughly in the range 
5-7). 

A similar analysis hold true for the other cases of small and large 612 , 613 , 
and thus will not be discussed further. To support the claim that the 
{P{i^e T^e)) surfaces for each value of $ are shifting in parameter space. 
Figure |^ is offered for comparison with Figure |^ This is the resulting surface 
for <I> = 4>sc- In addition, Figures present the case for = 0.8, s^^ = 
0.001. It is clear from both comparisons that the overall shapes of the sur- 
faces are the same, up to a shift of an overall constant. As suggested, this 
corresponds to an order of magnitude shift in the resonance boundary. 

Figures 11 and ^ are overhead projections of Figures |9|,|^ in the {2EAf2i, 
2£'A/3i) plane. The lighter-shade triangular-shaped wedge at the center 
depicts two of the features discussed here. The vertical boundary represents 
the resonance boundary, where {N^y^ > {Ne)max for each neutrino, while 
the diagonal one shows the intersection of the plane A/21 = A/31. It can be 
seen that the former shifts as discussed according to the potential while 
obviously the latter remains fixed. 



7 What's All This, Then? 

It thus becomes apparent that the same type of shifting effect as that noted 
in 0] could conceivably occur, depending on the choice of gravitational po- 
tentials affecting the VEP oscillations. Whereas before it was thought that 
VEP oscillations were singly dependent on the product $A/jj, these results 
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suggest that perhaps the value of the background potential <I> can indeed 
be constrained by the resultant neutrino data. The figures presented here 
immediately determine the observed energy spectrum of neutrinos, and 
many future neutrino observatories could detect these variations. As with 
the results in ||^, it is found that the shapes of the probabilities are 
strongly determined by the size of ^13, and small 13-angle oscillations pro- 
vide a much more diverse spectrum of {P{ve T^e)) curves for different 
input values (in this case, the potential <I>). 

Alternatively, what could also solve such a problem, should VEP be the 
mechanism at work, is the detection and subsequent spectrum analysis of 
extra-solar neutrinos from sources whose gravitational potential is known 
to either great accuracy or reliability. Various papers have addressed the 
detection of such high-energy intergalactic neutrinos subject to the VEP 



mechanism [21 



Additionally, the discovery of a possible "gravitationally- induced quan- 
tum phase" was recently discussed in the literature . While the reference 
treatment deals with MSW neutrinos which can possibly experience a phase 
shift due to interactions with strong gravitational sources, it may be possi- 



ble to extend the analogy to VEP neutrinos |23|. Should such an extension 
be possible, then it may be possible to determine the potential felt in the 
vicinity of the Sun, since the gravitationally- induced quantum mechanical 



phases of would be functions of the source potential (hence the product 
^fij^Source, and not just the product Afij^soiarSystem- These conclusions 
are as of yet unverified. 
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Figure Captions 

The following conventions are used for the figure conventions. The nota- 
tion df21, df31 on the plots corresponds to 2SA/2i,2£'A/3i which are 
expressed in units of 10^^ MeV~^. Also, P = {P{i'e '^e)), and s^^^^ = 
sin^ 26*12, 4i3 = sin^6'i3- 

Figure 1: (P(z/e ^ i^e)) curves for = 5 x lO'^sj^^ = 0.4, 2EAf3i = 

Figure 2: {P{ue ^ Ue)) curves for slg^^ = 0.8, s^^^ = 0.4, 2EAhi = 2x 10^, 
$ = <^0(r),00(i?0). 

Figure 3: {P{i^e ^ J^e)) curves for sjg^^ = 5 x lO'^, 'sj^^ = lO'^, 2EA/31 = 
2x 106, $ = 00(r),</.0(i?0). 

Figure 4: (P(i/e ^ i^e)) curves for = 0.8,5^^3 = IQ-^, 2£;A/3i = 

2x 106, $ = <^0(r),</.0(i?0). 

Figure 5: {P{ue ^ Ue)) curves for sjg^^ = 0.8, sj^^ = 10"^, 2£;A/3i = 
2 X 10^, $ = ^Q{r),^e{Re),(l)sc,(l>SD. 

Figure 6: (P(i/e ^ i^e)) curves for s^^^^ = 0.8,5^^3 = IQ-^, 2SA/31 = 
2 X 106, $ = ct>Q{r),(l>e{Re),(l>sc,(t>SD. 

Figure 7: (P(i^e ^ i^e)) surface for s^g^^ = 5x 10"^, Sg^^ = 10~^, $ = (pQ{r). 

Figure 8: (P(fe Ve)) surface for s^^^^ = 5 x 10"^, s^^^ = 10~^, $ = ^sc- 

Figure 9: {P{i^e i^e)) surface for s^g^^ = 5 x 10~^,Sgj3 = 0.4, $ = ^0(r). 

Figure 10: {P{i^e ^ z^e)) surface for s^g^^ = 5 x 10~^,s^j3 = 0.4, $ = ^sc- 

Figure 11: Overhead projection of Figure 9 showing resonance boundary 
and hierarchy breaking transition. 

Figure 12: Overhead projection of Figure 10 showing resonance boundary 
and hierarchy breaking transition. 
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